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Abstract. We show that intrinsic (not lensing-induced) correlations between galaxy shapes 
offer a new probe of primordial non-Gaussianity and inflationary physics which is complemen¬ 
tary to galaxy number counts. Specifically, intrinsic alignment correlations are sensitive to an 
anisotropic squeezed limit bispectrum of the primordial perturbations. Such a feature arises 
in solid inflation, as well as more broadly in the presence of light higher spin fields during 
inflation (as pointed out recently by Arkani-Hamed and Maldacena). We present a derivation 
of the all-sky two-point correlations of intrinsic shapes and number counts in the presence 
of non-Gaussianity with general angular dependence, and show that a quadrupolar (spin-2) 
anisotropy leads to the analog in galaxy shapes of the well-known scale-dependent bias in¬ 
duced in number counts by isotropic (spin-0) non-Gaussianity. Moreover, in the presence of 
non-zero anisotropic non-Gaussianity, the quadrupole of galaxy shapes becomes sensitive to 
far superhorizon modes. These effects come about because long-wavelength modes induce a 
local anisotropy in the initial power spectrum, with which galaxies will correlate. We forecast 
that future imaging surveys could provide constraints on the amplitude of anisotropic non- 
Gaussianity that are comparable to those from the Gosmic Microwave Background (CMB). 
These are complementary as they probe different physical scales. The constraints, however, 
depend on the sensitivity of galaxy shapes to the initial conditions which we only roughly 
estimate from observed tidal alignments. 
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1 Introduction 

The intrinsic alignments of red galaxies are thought to have their origin in the interaction of 
galaxy shapes with the tidal field of the large-scale structure [1] . These alignments remain a 
relatively unexplored cosmological probe, complementary to weak gravitational lensing and 
the statistics of the galaxy distribution. Upcoming imaging surveys, such as Euclid} and the 
Large Synoptic Survey Telescope (LSST^), will measure the shapes of billions of galaxies over 
half of the sky, providing a unique opportunity to exploit alignments as a cosmological probe. 

Alignments can contribute to the correlations between shapes and positions of galaxies 
with a strength larger than weak lensing at low redshifts, and sub-dominant at larger redshifts 
[2, 3]. Much of the effort in the study of alignments has been focused on measuring, modeling 
and removing intrinsic alignment contamination to weak gravitational lensing observables [4- 
7]. The tidal alignment model [1] has been shown to reproduce the scale-dependence of the 
observed alignment signal for red galaxies in the linear regime [8, 9] with an amplitude that 
depends on luminosity [10], while blue galaxies have no measured alignments and only upper 
limits for their intrinsic shape correlations are available [11, 12]. In this work, we use “red” 
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as a synonym of early-type and “blue”, of late-type galaxies. In the low-redshift Universe, 
a very small fraction of early-types have blue colors [13], but this assumption has not been 
tested at the high redshifts to be probed by upcoming imaging surveys. On the other hand, 
recent results from numerical hydrodynamical simulations suggest that indeed color is well 
correlated with galaxy dynamics as tracers of the intrinsic alignment signal [14]. Alignments 
are typically regarded as a contaminant to ‘cosmic shear’ (the two-point correlation of lensing- 
induced ellipticities), and they can contribute significantly to the correlation of galaxy shapes 
with the weak leasing of the Cosmic Microwave Background (CMB) [15, 16]. 

The exploration of alignments as a cosmological probe has only started recently. The 
wealth of shape information expected from upcoming imaging surveys, and the availability 
of a model that reproduces the observed alignment correlations of red galaxies, has triggered 
the exploration of the use of intrinsic alignments for cosmology. For example, the correlation 
between galaxy shapes and CMB B-mode polarization induced by primordial gravitational 
waves could be detected from future surveys, albeit with limitations in confirming its primor¬ 
dial origin [17]. While the CMB B-mode polarization is a purely linear effect induced by the 
gravitational redshift of tensor metric perturbations, the latter also induce an effective tidal 
field which modifies the local growth of structure and acts to align galaxies [18, 19]. Another 
example of the use of intrinsic alignments for cosmology is the possible detection of baryon 
acoustic oscillations in the cross-correlation of galaxy positions and shapes from future sur¬ 
veys [20]. While intrinsic alignments are weaker in amplitude than galaxy clustering, they 
do have some advantages such as the absence of redshift-space distortions (RSD) [21]. 

In this work, we show that intrinsic alignments can probe inflationary models where the 
squeezed-limit bispectrum of the primordial curvature perturbations, usually denoted TZ or 
is anisotropic. Phrased in terms of the primordial potential perturbation (j) during matter 
domination, which is related to the curvature perturbation by 7^ = (5/3)(/>, the squeezed-limit 
bispectrum can be expressed as 


B^(ki,k2,k3 




^ yl,n(ki • k5) P^{kL)P^{ks) 


1 + 0 



( 1 . 1 ) 


where = kL ki, k2 while k^ = ki — kx,/2, are the Legendre polynomials, and 
Ai are dimensionless amplitudes. Note that the angular dependence has to be even; see 
Section 2. In general, there can be several contributions with different power indices A. 
In this paper, we will focus on the “local” scaling with A = 0. However, our results are 
easily generalized to other scalings. In Eq. (1.1), the coefficient Aq of the isotropic term is 
related to the usual local non-Gaussianity parameter via Aq = 4/^£. We are interested 
in constraining the parameter that governs the leading anisotropic contribution, namely the 
quadrupolar dependence of the bispectrum, A 2 . As we will show, A 2 does not lead to a 
significant large-scale effect in galaxy clustering, specifically the two-point function. It does 
however lead to the analog of the well-known scale-dependent bias for galaxy counts, which 
scales as f^^{k/aH)~'^, in the two-point statistics of galaxy shapes, adding a contribution oc 
A 2 {k/aH)~‘^. Combining galaxy clustering with shapes, we can then constrain both isotropic 
and anisotropic non-Gaussianity simultaneously. 

Thus, alignments offer a golden opportunity to explore a new dimension in the parameter 
space of inflationary models. Anisotropic non-Gaussianity can arise in several early-Universe 
scenarios. In solid inflation [22], inflation is driven by the exponential stretching of an 
unusual solid. A feature of solid inflation is that it produces large anisotropic non-Gaussianity 
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while is small {A 2 S> ^ 0 ) [22]. Both isotropic and anisotropic non-Gaussianity are 
produced by primordial curvature perturbations generated by large-scale magnetic fields 
[23, 24]. Anisotropic non-Gaussianity can also be produced in inflationary models with a 
generalized bispectrum from excited Bunch-Davies vacuum [25]. Recently, Arkani-Hamed 
and Maldacena [26] pointed out (see also [27-29]) that the squeezed-limit bispectrum of 
primordial perturbations offers a clean probe of the spectrum of particles with masses of 
order the Hubble scale and smaller during inflation (see [30] for a corresponding study of 
the tensor-scalar-scalar three-point function). The index A is related to the mass of the 
particle, with A —)• 0 for degrees of freedom much lighter than the Hubble scale. The order 
of the angular dependence L encodes the spin of the particles, as the leading contribution of 
a spin-s particle yields a contribution of order L = s in Eq. (1.1). Thus, our results show 
that the combination of galaxy clustering and shapes can probe the presence of scalar (spin 
0) and tensor (spin 2) degrees of freedom with m < H during inflation. 

Of course, this type of primordial non-Gaussianity can also be constrained using the 
CMB, which yields the currently tightest constraints of A 2 = 13 ± 93 (la) [31]. We stress, 
however, that galaxy clustering and shapes probe the three-point function in Eq. (1.1) for 
small-scale modes that are on significantly smaller scales than the GMB. Given that the 
amplitude of non-Gaussianity can well be scale-dependent, the constraints we forecast are to 
be seen as complementary to the CMB. 

Our work is organized as follows. Section 2 presents the non-Gaussian matter den¬ 
sity field and our conventions for parametrizing anisotropic non-Gaussianity through the 
squeezed-limit bispectrum. The impact of tidal fields and primordial non-Gaussianity on 
the shapes and clustering of galaxies is described in Section 3. In Section 4, we derive full 
sky observables of intrinsic alignments. Section 5 summarizes the adopted model for the red 
galaxy population, based on [3], while we outline the Fisher forecast in Section 6. In Section 
7, we present a forecast of the detectability of anisotropic non-Gaussianity in an LSST-like 
weak leasing survey. Finally, we discuss our results in Section 8. 

Throughout this paper, we assume the following Planck [32] fiducial flat ACDM cos¬ 
mology: = 0.022, = 0.12, h = 0.67, Q.k = 0, Ms = 2.2 x 10~®, Ug = 0.9645, 

kp = 0.05 Mpc“^ and we define = H;, -|- Hcdm- 

2 Large-scale non-Gaussianity 

When computing two-point correlations of large-scale structure in the large-scale limit, the 
leading relevant non-Gaussian statistic is the squeezed-limit potential bispectrum ^^(ki, k 2 , k^), 
where kr ^ ki, k 2 . While bispectra are in general non-separable functions of the momenta, 
they can usually be well approximated by a sum of separable terms [33]. For simplicity, we 
will only consider one of those terms, which in the squeezed limit can be written in general 
as 

H^ki, ks, ks) = f{kL, ks)P^{kL)P^{ks) , (2.1) 

with “beyond-squeezed limit” terms suppressed by (ki/ks)'^. These are negligible when 
considering large-scale correlations, and we will not consider them here. If the bispectrum is 
scale-invariant, which we will assume in the following, then f{kL,ks) = /(k^ • k 5 , k^/ks). 
Further, one can parametrize the dependence on ki/ks as a power law (ki/ks)'^, so that 
/(k/,,k 5 ) —)• /(ki • ks) (kL/ks)^. In the following sections, we will restrict to bispectra of 
the local-type for simplicity (A = 0). The generalization beyond the simplifying assumptions 
made here is straightforward and will be pointed out below where relevant. 
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Finally we can decompose the angular dependence /(/u) into multipole moments, where 
// = ki • ks: 

f{l^)= ( 2 - 2 ) 

£= 0 , 2 , 4 ,--- 

It will become clear below why we neglect odd multipole moments. This then leads to 
Eq. (1.1). Standard local non-Gaussianity corresponds to L = 0 and Aq = 4/^l- We are 
interested in the two lowest order contributions: a monopole Aq, and a quadrupole A 2 , which 
corresponds to the anisotropic shape generated, for example, by solid inflation [34] or spin-2 
particles [26]. Thus, /^l = 0 or “no local non-Gaussianity” will correspond to the statement 
that ^0 = 0 in Eq. (2.2), while A 2 = 0 corresponds to a vanishing quadrupole. 

To interpret the effect of the bispectrum defined in Eq. (2.1), consider a single long- 
wavelength plane wave potential perturbation 4>{kL). According to Eq. (2.1), the local small- 
scale power spectrum at position x is given by 


(k5,x|())(kL)) 


1 + E 


£= 0 , 2 , 4 ,- 




P<p{ks) 


(2.3) 


where P{ks) is the Gaussian power spectrum. This can be easily verified by multiplying 
by <f>{k'j^), integrating over d^x, and taking the expectation value. Since (/> is a real field, 
P(j){—ks) = P^{ks) has to hold, and Eq. (2.3) implies that only even multipoles can contribute 
to the leading squeezed-limit bispectrum Eq. (2.1). 


3 Intrinsic alignments and non-Gaussianity 

In this section, we derive the linear relation between intrinsic galaxy shapes and the under¬ 
lying matter distribution in the presence of primordial non-Gaussianity, including the effect 
of the anisotropic bispectrum. Intrinsic alignments can be treated in exactly the same way as 
galaxy biasing, the only difference is that projected shapes are a spin-2 quantity on the sky, 
as opposed to the scalar galaxy density. Thus, our discussion parallels that of [35-38] who 
studied the impact of primordial non-Gaussianity on galaxy clustering. We describe galaxy 
shape statistics in two stages. First, we provide a description of the statistics of the physical, 
intrinsic three-dimensional shapes of galaxies, denoted as gij. We will assume that gij is 
related to the trace-free part of the second moment tensor lij of, for example, the intrinsic 
emissivity of a given galaxy; 

so that gij is dimensionless. Second, we project the shapes onto the sky plane, thus providing 
the connection to the observed two-dimensional shapes, and add the contribution due to 
gravitational lensing. Since the processes that determine galaxy shapes are intrinsically 
three-dimensional, we argue that a physical treatment of alignments should proceed in this 
way, rather than starting with an expansion on the projected galaxy shapes. 

The precise location and shape of a given galaxy is determined by the nonlinear distribu¬ 
tion of matter within some finite spatial region around the galaxy (along its past trajectory) 
described by a typical scale i?*. For example, if the dominant mechanism for determining 
position and shape is gravitational collapse, i?* will be of the order of the Lagrangian ra¬ 
dius of the parent halo of the galaxy. Since we are interested in the large-scale statistics of 


9ij — [44 -^fci] 


-1 


hj 2 ^ij 


(3.1) 
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galaxies, let us coarse-grain the galaxy distribution on some scale much larger than (but 
smaller than the scale of correlations we are interested in). Then, the only properties of the 
matter distribution that are relevant to galaxy formation are the matter density 5 = p/p — 1 
and tidal field Kij, coarse-grained over a scale i?* (e.g., [39, 40]).^ We define the tidal field 
via 

^ dvrG p a? 


didj - 


^ = Vij5, 


(3.2) 


where is the gravitational potential and 


A. = ^ • (3.3) 

The quantity Kij is trace-free, so that the six degrees of freedom contained in {h, Kij} fully 
describe the second derivative tensor didj^. Due to the equivalence principle, an observable 
such as galaxy shapes cannot depend on or itself. 

Thus, on large scales, r S> i?*, the statistics of galaxy shapes and number, quantified 
by the fractional overdensity 6n, can be effectively described through local functions of the 
coarse-grained 6 and Kij: 


h„(x,T) =T[(5(x,T),K'fci(x,r),r] 

9ij{^,T) = Gij[h(x,r),Kfcz(x,r),r], (3.4) 


where r denotes conformal time. The leading correction to Eq. (3.4) due to the non-locality 
of galaxy formation is given by terms of order and which we will neglect 

throughout this work. In Fourier space, they lead to a scale-dependent bias oc k'^R^ which 
becomes relevant on small scales. 

In Eq. (3.4), we have ignored the small-scale perturbations which are also present in 
the real universe. This is fine in the case of Gaussian initial conditions, since the small-scale 
perturbations are not correlated with the long-wavelength perturbations [apart from terms 
that already appear in Eq. (3.4)]. As we will see, this will not be the case for non-Gaussian 
initial conditions. 

A Taylor expansion of Eq. (3.4) then leads to the usual local bias expansion for the 
galaxy number density perturbation 6n, i-e. 

h„(x,r) = 6)*(r)(5(x,r) + ^^^(r) ^^(x, r) + ^6 )^(t) (Kijf(^,T) + ••• . (3.5) 

The analogous expansion for pij, up to quadratic order, leads to 




gij{x,T) =b{{T)Kij{x,T) -F -b 2 {T)Kij{x,T) 6 {x,T) + -btir) 


KikK) - MKirn? 


+ ■■■ 


(x, t) 
(3.6) 


Note that we have allowed for all possible terms consistent with the symmetries of 6n and 
Pij, and in the absence of preferred directions in the galaxy frame. The coefficient 6” cor¬ 
responds to the z-th order response of the mean number density of galaxies to a change in 

^More accurately, the relevant properties are the coarse-grained density and tidal field along the entire 
past trajectory of the fluid elements. However, this is not relevant at the order up to which we work here. 
See [41, 42] for a discussion in the context of galaxy bias. 
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the background density, i.e. bf = {p^jng)d'^rigidp^. Similarly, b{ corresponds to the linear 
response of the mean shape of galaxies to an external tidal field Kij. Due to a lack of pre¬ 
ferred directions in the frame of the galaxy, this response must have the simple form given 
in Eq. (3.6). That is, b[ is simply a number rather than a tensor. Similarly, b^ corresponds 
to the response of galaxy shapes to an external tidal field and a simultaneous change in 
the background density. Such a term is expected to be present intrinsically, however it is 
also induced by the weighting by the galaxy number density [43]. Finally, b^ quantifies the 
quadratic tidal alignment. We emphasise that Eqs. (3.5)-(3.6) are effective relations describ¬ 
ing galaxy statistics on large scales rather than deterministic relations valid at a given scale. 
We have also neglected stochastic terms in the bias relations, which have to be included in 
general. The cubic and higher order terms neglected in Eqs. (3.5)-(3.6), while important on 
small scales, are not relevant for the discussion here as they do not lead to scale-dependent 
signatures of primordial non-Gaussianity on large scales. 

Let us now consider the cross-correlation of galaxy shapes with the matter density held 
on large scales. Eq. (3.6) yields 


{S{yi)gij{y)) = 6f(<5(x)Ky (y)) -bi {6{^)6{y)Kij{y)) 


+ 2^1 i-r) ( <^(x) 


KikK’^i - -SiffKirn? (y) 


(3.7) 


where the ellipsis stands for higher order correlators which are not relevant on the large scales 
we are interested in. The hrst term in Eq. (3.7) is the leading term in the case of Gaussian 
initial conditions. It is given, in terms of the matter correlation function ^(r), by 


{6{:x.)Kij{y)) = AjCdx - y|). (3.8) 

The second and third terms in Eq. (3.7) involving three-point functions are exclusively in¬ 
duced by non-linear evolution in the case of Gaussian initial conditions, in which case they 
are only relevant on small scales. Hence, the large-scale limit of the matter-shape cross¬ 
correlation becomes 

(5(x)g,,(y)) - y|) , (3.9) 

or equivalently in Fourier space, 

((5(k)ffij(k')) = b{ Pynik) (27r)3(5£,(k -h k'), (3.10) 

where Pm{k) is the matter power spectrum. This implies that b{ is the analog of the usual 
linear bias for galaxy density correlations. Eq. (3.10) corresponds to the prediction of the 
so-called linear alignment model [Ij. This is generically the leading contribution to galaxy 
shape correlations on large scales, unless b{ happens to be very small for a given galaxy 
sample. 

We now consider the case of non-Gaussian initial conditions of the local type, quantihed 
by Eq. (1.1). The three-point terms in Eq. (3.7) are sensitive to the bispectrum in the 
squeezed limit, since two fields are evaluated at the same point. In the following, for simplicity 
of the presentation we will only focus on the first term oc b^', the term oc bj- behaves in 
essentially the same way (see App. A), so including it here would add nothing new. As 
shown in App. A, for a squeezed-limit bispectrum of the form Eq. (1.1) we have 

{6{x)6{y)Kij{y)) = - y\){S^), (3.11) 
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where ^ 5 (f,{r) is the linear matter-primordial potential cross-correlation. Equivalently, in 
Fourier space, 

(,5(k)[5Ki,](kO> = ^A 2 M-\k)Pm(k)(d^} (27rfdn(k + k'), (3.12) 

where Pm{k) is the matter power spectrum in synchronous-comoving gauge, and 

_ 2pT{k)D{z) 

~ 3 

is the linear relation between primordial potential and density, (5(k, z) = M.{k, z)(p{k), with 
T{k) denoting the transfer function and D{z) denoting the growth factor normalized to 
D{z) = 1/(1 -|- z) during matter domination. Thus, there is a non-Gaussian contribution to 
the large-scale shape correlation if the squeezed-limit bispeetrum is anisotropic, A 2 7 ^ 0 . 

However, Eqs. (3.11)-(3.12) involve the variance of the matter density held (<5^), which is 
formally divergent if we let the smoothing scale go to zero, and hence needs to be regularized 
by introducing a counter-term. Equivalently, it corresponds to an unphysical dependence of 
our observable on the arbitrary smoothing scale. We see that the counter-term we need to add 
to the bias expansion of galaxy shapes, Eq. (3.6), is Vijcp. This term is not due to gravitational 
evolution, since it involves the potential itself and does not obey the equivalence principle, 
but is due to the initial conditions. Note that in Eq. (3.4) we have implicitly assumed that 
the statistics of the initial small-scale modes which inhuence local galaxy formation are the 
same everywhere. The squeezed-limit bispectrum of Eq. (1.1), however, precisely describes 
a modulation of the power spectrum of small-scale density fluctuations by long-wavelength 
potential perturbations. This exactly parallels the case of galaxy density correlations in the 
case of local-type primordial non-Gaussianity [35]. The key difference for the generalized 
squeezed limit of Eq. (1.1) is that now the local power spectrum becomes anisotropic. Thus, 
in the case of primordial non-Gaussianity with ^2 7 ^ 0 we have to explicitly take into account 
the dependence of galaxy shapes on the amplitude and anisotropy of small-scale fluctuations. 

Consider the following transformation of the initial matter power spectrum, 

= (1 + 2ai^Uk^) P^lsoik) (3.14) 



(3.13) 


where aim is trace-free and Pm\so i® fiducial isotropic power spectrum. We can define 
the mean response of galaxy shapes to such an anisotropic initial power spectrum. 


^NG 


9{9ij)a 

doij 


Q ;=0 


(3.15) 


where no summation is to be taken (the absence of other preferred directions requires b^Q to 
be a number rather than a tensor). As we show in App. B, allowing for this dependence of 
galaxy shapes on the anisotropy of the initial power spectrum in Eq. (3.4) provides the correct 
counter-term 3b^QA2'Dij(j) to cancel the divergence appearing in Eqs. (3.11)-(3.12) (the same 
holds for the divergence appearing in the other three-point term oc bj-). This is not surprising, 
since according to Eq. (2.3) (with A = 0) a long-wavelength potential perturbation yields a 
modulation of the type Eq. (3.14) with 


2a'™ = A 2 




4>(kL)P^^^ = A2P'™</.l(x) . 


(3.16) 
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The end result is that the large-scale matter-shape cross-correlation becomes 

{6{x.)gij{y)) = - y|) -h - y|), (3.17) 

and equivalently in Fourier space 

{6{k')gij{k)) = l^b{ + 3bi^QA2 Pm{k) {27rf5D{k+ 'k'). (3.18) 

This result is very similar to the well-known scale-dependent bias induced by local 
primordial non-Gaussianity in the two-point function of the galaxy density. Specifically, the 
galaxy density-matter cross-correlation is given by [37] 

(5(k')<5„(k)) = !^b^ + ^b'^(.AoM-\k)^Pmik)i27Tf6Dik + k'), (3.19) 

where b^Q = n~^dfig/d In ag is defined as the response of the galaxy abundance to a change 
in the (isotropic) initial power spectrum amplitude quantified by ag. For dark matter halos 
and assuming a universal mass function, we can relate this to the linear density bias via 
&NG = 

Two facts are noteworthy about Eqs. (3.18)-(3.19). First, isotropic local type non- 
Gaussianity {Aq / 0; ^2 = 0) does not generate any large-scale correlations of galaxy shapes 
at linear order in the non-Gaussianity. The reason is that, due to symmetry, an isotropic 
rescaling of the small-scale fluctuations cannot affect galaxy shapes. Second, there is no 
impact of anisotropic non-Gaussianity (Aq = 0; ^2 / 0) on large-scale galaxy density cor¬ 
relations. The number density of tracers depends on the amplitude of small scale fluctua¬ 
tions. It can be shown (App. A) that this is not modified by anisotropic non-Gaussianity. 
Thus, at the level of two-point functions, only galaxy shapes are affected by anisotropic 
non-Gaussianity, while the same is true for the galaxy density and isotropic non-Gaussianity. 
By comparing shape correlations with density correlations we can thus disentangle isotropic 
from anisotropic shapes of the local type. 

Finally, it is straightforward to generalize Eq. (3.19) to different squeezed-limit power 
scalings A. Specifically, 

(,5(k0ff.,(k))= + >f-i(fc)|p^(fc)(27r)3<5^(k + k'), 

(3.20) 

where kp is a pivot scale and we have indicated that now ^ depends on A (and kp). Eor 
A ~ 2, the non-Gaussian contribution is only weakly scale-dependent and thus becomes more 
difficult to constrain. 

To summarize, the large-scale correlation between matter and galaxy shapes has a Gaus¬ 
sian piece oc b {, where b{ represent the linear response of galaxy shapes to a long-wavelength 
tidal field; and a non-Gaussian piece, which is proportional to the quadrupolar part of the 
squeezed-limit bispectrum. The coefficient of the non-Gaussian term, b^Q, quantifies the 
linear response of galaxy shapes to an anisotropic power spectrum of initial fluctuations 
[Eq. (3.14)]. Our derivations in this section represent the generalization to a spin-2 held of 
the well-known non-Gaussian halo bias [35, 37, 44]. 

While it is clear that the scale-dependent bias of galaxy shapes has to be present in 
general, it is difficult to predict a precise number of b^Q for dark matter halos, let alone 
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Figure 1. The tidal bias induced by anisotropic non-Gaussianity, Eq. (3.21), for b{ = 
-0.in^D{0)/Diz), = -0 • irim and A 2 = —10 (see Section 5). The increasing thickness of 

the curves represents different redshifts from z = 0 to z = 1.5. At fc > 10“^ h Mpc“^, the effect of 
non-Gaussianity is negligible and all curves asymptote to —6{. Anisotropic non-Gaussianity changes 
the tidal bias on large scales and the effect is larger at high redshift at a given scale. 

galaxies. If it is very small, then this effect becomes less interesting observationally. The 
value of b^Q could be measured by running N-body simulations with anisotropic initial power 
spectrum and measuring the mean shapes of halos along the anisotropy axes. This is beyond 
the scope of this paper. In the context of the excursion set [45] however, one would expect 
that b^Q is comparable to the linear tidal alignment coefficient b{. In both cases the collapse 
threshold u = becomes direction-dependent, i.e. different along different axes. In the 
tidal case, this is due to the faster infall along the steeper potential gradient (e.g., [46]), while 
an anisotropic power spectrum of small-scale fluctuations will lead to a similar anisotropic 
V. For our forecasts, we will thus assume that ~ b{. This should only be understood as 
an estimate at the order-of-magnitude level however. That is, while we expect the strongest 
signature of anisotropic primordial non-Gaussianity to appear for those galaxies which show 
the strongest linear tidal alignment (i.e. massive elliptical galaxies), the amplitude of the 
two coefficients could certainly differ by a factor of several. 

Figure 1 shows the total effective bias of intrinsic alignments, 

b{+ 3bf^QA2M~^{k,z), (3.21) 

i.e. the term in curly brackets multiplying Pm{k) in Eq. (3.18), as a function of scale. In 
the Gaussian case, it is a constant, while for A 2 7 ^ 0 it is strongly scale-dependent. Here, 
we assume a fiducial value of b{ = (consistent with LRG measurements 

at low redshift [ 8 ]), and = —O.lHm (see Section 5). The impact of non-Gaussianity 
increases with redshift. Depending on the sign of A 2 , the anisotropic tidal bias can have the 
same sign as the isotropic linear alignment model (if A 2 > 0 ) or the opposite sign (if A 2 < 0 ). 
In the latter case, the effect of the anisotropic tidal bias is to generate tangential alignments 
on large scales. This transition scale can be estimated by simply setting Eq. (3.21) to zero, 
and it is A: ~ 1.5 x 10“^ h/Mpc for A 2 = —10. The impact of this bias on angular shape 
correlations will be shown in the next section. 
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4 Observed shape and number count correlations 

In this section we derive the implications of the scale-dependent bias in galaxy alignments for 
observed angular correlations of galaxy shapes and number counts. We will also include the 
contributions from gravitational leasing to shapes, which are comparable to or larger than 
alignments at cosmological redshifts. However, we will neglect the impact of primordial non- 
Gaussianity on the leasing correlations, since they are highly suppressed by the projection 
along the lightcone [47]. This is discussed in detail in App. C. 

4.1 Projection 

In imaging surveys, one does not have access to the 3D shape of galaxies which we have 
written in Eq. (3.17). Instead, one only measures (in the absence of leasing) the intrinsic 
shape projected onto the sky: 


7/,*i 





kl 


9kl , 


(4.1) 


where 

T^ij ~ ^ij fliflj (4.2) 

is the projection operator and h is the unit vector along the line of sight. Note that VijV^^ = 
pA, and any term oc Ski in 9ki drops out in the projection Eq. (4.1). In Eq. (4.1) we have 
introduced the notation for the intrinsic shape, as we will also encounter the leasing 
contribution 7G,ij below. The observed shape is the sum of both contributions. 

Further, in imaging surveys galaxies follow a redshift distribution, dN/dz. Thus, fol¬ 
lowing Eq. (3.17), the intrinsic alignment contribution to observed shapes of galaxies in the 
direction h are related to the density field and gravitational potential through 




(h) = j dz 


dN 

dz 


t{z)) + 2>bl^QA2Vl^(j){x{z)rv] t{z)) 


(4.3) 


where xi^) is the comoving distance out to redshift z, t{z) is the conformal time at that 
redshift, and we have defined the projection onto the sky of the derivative operator Dij 
through 

^ = (Vi Vki . (4.4) 

Given the discussion of the renormalization approach in Section 3, Eq. (4.3) is not to be seen 
as a truly local relation, but merely as formal relation in order to derive angular correlations. 

Eq. (3.18) shows that the effects of primordial non-Gaussianity on alignment correlations 
are relevant on large scales. For this reason, we calculate proper full-sky angular correlations. 
We begin by reviewing the formalism employed in the following to obtain full-sky correlations 
efficiently. 7 ij(h) is a traceless 2-tensor on the sky, with two degrees of freedom. Following 
the literature developed for GMB polarization [48-50], we decompose Xij components 
that transform with a spin s = ±2. Gonsider an orthonormal coordinate system (ei,e 2 ,n). 
If we rotate the coordinate system around h by an angle ij}, so that —>• e(, the linear 
combinations = (ei =p ie 2)/\/2 transform as 


m± —)• m(j_ = 


(4.5) 
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A general function /(n) is spin-s if it transforms under the same transformation as 


fih) ^ /(fi)' = e“^/(A). (4.6) 

An ordinary scalar function on the sphere is clearly spin-0, while the unit vectors defined 
above are spin±l fields. This decomposition is particularly useful for deriving multipole 
coefficients and angular power spectra. For a summary of useful results, see App. A of [51]. 
We also define 


k± = 

E± = (4.7) 

for any 3-vector fcj and 3-tensor Eij. We then decompose (this equally applies to the 
individual I and G contributions) as 

lij = 27 + _27 

±27 = , (4.8) 

where ±27 are spin ±2 functions on the sphere (in analogy to Stokes parameters Q ± ill). 
Consider the contribution of a single Fourier mode k of the density and potential, which we 
choose along the 2 axis, to ± 27 /(h): 

± 27 /(h,/cz) = J dz^ (x) + 3^NG^2<^(k;T(2:))] 

= - ^(1 - /^^) J dz^ [bidmO^; t{z)) + 36^G^2<(>(k; t{z))] , (4.9) 

where /r is the cosine between k and h, and we have defined x = kx{z). We have also used 
that Vijm^^rrij. = 0. Since ±27 are not scalar quantities on the sky, they do not transform 
trivially under rotations. Instead, we apply spin-raising 9^ and lowering 9^ operators to 
obtain a scalar on the sky (see App. A of [51] for details). Specifically, for a spin -2 quantity 
such as ± 27 ; the spin-lowering operator is given by [see Eq. (A 8 ) of [51] for m = 0] 

Q 2 

8 ^ ±27/ = ^ [(1 “ //^)±27/(n,z)] 

= ^ / dz^ [6{5m(k; t{z)) + 36^GA2(/>(k; r(z))] (1 5^)^ [x^e^^] . (4.10) 

In the last line, we have used 

^ [(1 - ±)±“'‘] =1^ [(1 + aS±”'‘] = -(1 + 83= . (4.11) 

This derivative operator is real, which reflects the fact that the linear alignments from scalar 
perturbations considered here do not produce S-mode shape correlations even in the presence 
of anisotropic non-Gaussianity. 5-modes will be generated on small scales by the nonlinear 
alignment terms written in Eq. (3.6), however. Applying the spin-raising operator 9^ to - 27 / 
yields exactly the same result. Crucially, 9^ ± 27 / as well as 9 ^_ 27 / transform as scalars on 
the sky, which allows us to evaluate their two-point function (angular power spectrum) in a 
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straightforward way. The angular power spectrum of 9^ + 27 / is then directly proportional to 
that of 7 /, which is our final goal (see Section 4.2). 

We emphasize that Eq. (4.10) is the exact, proper relativistic expression of intrinsic 
shapes on linear scales. This is because, as shown in [40], the tidal field defined in conformal- 
Newtonian gauge corresponds to the physical, locally observable tidal field at linear order. 
The leasing contribution discussed in Section 4.4 below then corresponds to the mapping 
of the intrinsic locally observable shape to our observations from Earth. Note also that 
redshift-space distortions do not affect alignments at linear order, since at that order they 
only modify the physical volume which does not affect shapes (see also [21]). 

Correspondingly, we can also write down the observed galaxy overdensity at linear order: 


5n(h, kz) = j 


b^6m{k;T{z)) + -b^QAo<p{k;T{z)) 


^IX/l 


(4.12) 


While the expression for the intrinsic shapes Eq. (4.10) is exact in the linear regime, we have 
neglected several linear-order contributions in the galaxy overdensity: 


• Redshift-space distortions: We will consider photometric galaxy samples throughout, 
for which RSD are a subdominant effect. Eurther, at / = 2 it is inconsistent to only 
include the linear RSD term without further relativistic corrections (see below). 

• Magnification bias: This contribution, while straightforward to include, is typically sub¬ 
dominant and furthermore affected negligibly by primordial non-Gaussianity as shown 
in [47] (see also App. C). 

• Relativistic corrections: There are corrections to the galaxy clustering kernels (e.g., 
[52] and App. B in [53]) that scale as (aH/k) and {aH/kfi. The latter terms are of the 
same form as local-type non-Gaussianity with Aq < 1 for typical sample parameters. 


These simplifications are justifiable, since the focus of this paper is shape correlations rather 
than clustering, and since these approximations are not expected to have a signihcant impact 
on the forecasted constraints for A 2 . 

4.2 Intrinsic shape correlations 

We can now apply the results of App. Al of [51] to obtain the angular power spectra of galaxy 
overdensity and shapes, and their cross-spectrum. In the following, we will consider three 
different redshift distributions of galaxies in order to keep the treatment general: dNc/dz 
denotes the redshift distribution of all galaxies with measured shapes; dN^-e^/dz denotes 
the subset of these galaxies that are early-type and known to show linear alignment of the 
type considered here (a linear alignment has so far not been detected for late-type galaxies). 
Einally, dN^/dz denotes the galaxy sample used in the measurement of the clustering. Eor 
our main results, we will assume that Nn, A^red and Nq are the same population of red 
galaxies. We summarize our model for this population in Section 5. 

Eq. (4.10) is a special case, r = 0, of Eq. (A19) in [51], and Eq. (A24) there (for |s| = 2 
and A^p = 1, where Np denotes the number of polarization states) immediately yields the 
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(a) A 2 = 10. (b) A 2 = -10. 


Figure 2. The angular auto-power spectrum of intrinsic shapes (without accounting for weak grav¬ 
itational leasing) for our fiducial choice of and A 2 = -1-10 (left) and for A 2 = —10 (right). The 
Gaussian case is shown in solid black for three redshifts, 2 = {0.5,1,1.5} from top to bottom. We 
show the case with anisotropic non-Gaussianity at the same redshifts (dashed). 


auto-power spectrum of the intrinsic component of the shapes: 


C'/(0 = ^|q:|y / k'dk^rn(k)lF/(k)l^ 


Fl{k) = [b{ + 3b^^QA2M ^{k, z)] [(1 + S^)^ x'^ji{x)]^^ 


\ x=kx(z) 


(^ + 2) ! 

(1-2) 


! / [6( + 36^g^2-M \k,z)] 


ji{x) 


x=kx(z) 


(4.13) 


(4.14) 


Here, Pmik) is the matter power spectrum at z = 0, and M.{k,z) is defined in Eq. (3.13). 
ji denotes the spherical Bessel function of order 1. Note that no B-modes are generated 
by linear scalar contributions. Also, by symmetry the lowest non-zero multipole of shape 
correlations is the quadrupole 1 = 2. 

Similarly, we can write down the cross-correlation between galaxy overdensity and shape 
E-modes (again, only the intrinsic part). This is the multipole space version of the density- 
tangential shape correlation function in real space, and reads 


Cnl{l) = - 

TT 

FPik) = J dz 


2 /(^-2)! / ,2 


ttV (^ + 2)! J 

dNn D{z 


k^dk P„,{k)Fl {k)Fp{k) 


dz D{0) 


bi + -6ng^oA4 {k, z) 


ji{x) 


x=kx{z) 


The clustering auto-power spectrum, on the other hand, is simply given by 


Cnn{l) = - / k^dk Pm{k)\FP{k)\^ . 
TT 7 


Eqs. (4.13)-(4.17) are valid on the full sky. 


(4.15) 

(4.16) 


(4.17) 
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4.3 Signatures of non-Gaussianity 

We now illustrate the effect of isotropic and anisotropic non-Gaussianity on the angular power 
spectra of galaxy counts and the intrinsic contribution to shapes (the leasing contribution 
will be considered in the next section). Depending on the sign of and A 2 , we obtain an 
enhancement or suppression of the angular correlations at low multipoles (as expected from 
Figure 1). In this section, we will assume the same values for the bias parameters as for 
Figure 1, i.e. b( = —0.1Qm^(0)/D(z) and = —O.lDm- 

By definition, the intrinsic shape auto-power spectrum (Figure 2) is positive regardless 
of anisotropic non-Gaussianity. This power spectrum decreases with redshift, as expected 
from Eq. (4.20), from 2 : = 0.5 (top) to 2 ; = 1.5 (bottom). Anisotropic non-Gaussianity 
changes the power spectrum for small multipoles (I < 20) and asymptotes to the Gaussian 
case on small angular scales. The fractional effect increases towards higher redshifts. 

An interesting and novel effect of primordial non-Gaussianity appears in the quadrupole 
of shape correlations for A 2 7 ^ 0: consider the term oc (^ 2 )^ in Cj{l). Focusing just on the 
/c-dependence at an effective source redshift Zg, this contribution to Cj{l) is proportional to 


Ci{l) oc j k^dkP^{k){bl^GA 2 


ji{x) 


x=kx{zs) 


(4.18) 


Now consider the quadrupole I = 2. We see that the factor in brackets asymptotes to a 
constant for kxizg) ^ 1 , and this term becomes, for k <C l/x{zg), 


Ci{l) {bi,GA2f J k^dkP^k). (4.19) 

This quantity diverges logarithmically in the IR for a near scale-invariant spectrum of primor¬ 
dial fluctuations. Thus, the quadrupole of galaxy shapes depends on the lowest wavenumber 
^min at which non-zero anisotropic primordial non-Gaussianity is present, no matter how 
small. In principle, a measurement of Cj{l = 2), together with the higher I which con¬ 
strain A 2 , would thus allow us to measure this low k cut-off, which is presumably related 
to the beginning of inflation (fcmin ~ 1 /Tinf where rinf is the conformal time at the begin¬ 
ning of inflation). Quantitatively, the contribution of such super-survey modes is roughly 
of order {A 2 )‘^As\'n[Ti-a_i/x{zs)]-, where we have neglected the contribution from the small tilt 
Ug — 1. Given the smallness of the amplitude of fluctuations Ag ~ 10“^^, this is probably 
only observable for an extremely long period of inflation with correspondingly large r^nf (or, 
alternatively, if the fluctuations near the begining of inflation happened to be much larger 
than the extrapolation of the near scale-invariant spectrum). For definiteness, we choose 
fc min = 2 X 10“®/iMpc“^ for our numerical results; for A 2 = 10, Cj{l = 2) decreases by 
roughly a factor of 3 when increasing fcmin to 10“^ /iMpc“^. 

It might be surprising that modes that are far outside our current horizon affect observ¬ 
ables. Physically, this is because for A 2 A 0) the matter power spectrum in our observable 
universe is anisotropic and has a spatially constant quadrupole with an expectation value 
proportional to (^ 2 )^/ k‘^dkP^{k). Galaxy shapes respond to this anisotropy, as quantified 
by b^Q, and the spatially constant power spectrum quadrupole will lead to a correspond¬ 
ing quadrupole of galaxy shapes. Related effects of superhorizon modes have been found in 
[54, 55]. Note that a similar effect also exists for isotropic non-Gaussianity Aq A 0, where the 
monopole of the galaxy overdensity similarly receives a contribution oc (Aq)^ f k‘^dkP^{k). 
However, this effect is unobservable, since the amplitude of scalar fluctuations themselves 
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(a) Positive A 2 combined with Aq = {0, 10} and com- (b) Negative A 2 combined with Aq = {0,10} and 
pared to the Gaussian case. compared to the Gaussian case. 


Figure 3. Fractional change in the angular cross-spectrum of galaxy positions and intrinsic shapes 
at z = 1 for different levels of non-Gaussianity and our fiducial choice of b^Q. We assume &" = 2 and 
^NG = (^1 ~ The left panel shows results for A 2 = -|-10, while the right panel shows A 2 = —10. 
In addition, we show results for isotropic non-Gaussianity, A 2 = 0, Aq = -1-10 in both panels. 


receives the same correction, which is thus absorbed in the normalization As- In this regard, 
anisotropic PNG and galaxy shapes in principle offer qualitatively new information over the 
well-known isotropic case. 

When the cross-spectrum of galaxy positions and intrinsic shear is considered, both 
A 2 and ^0 have an impact on the large-scale correlations. Figure 3 shows the effect of the 
different types of non-Gaussianity on this cross-spectrum and assuming 6” = 2. For A 2 < 0 
(right panel), the correlation between galaxy positions and shapes changes sign on large 
scales, corresponding to a transition from radial to tangential alignment. That is, intrinsic 
shapes in this case have the same sign as the lensing signal, which is also tangential behind 
overdensities. If A 2 > 0, non-Gaussianity strengthens the effect of the tidal field, enhancing 
the radial alignment of red galaxies. In this case, there is no sign change. We emphasize that 
both A 2 and ^0 are theoretically allowed to have either positive or negative signs. 


4.4 Lensing shear 

Gorrelations in galaxy shapes are not only present intrinsically through tidal alignments, 
but are also induced by gravitational lensing (in particular, the shear; see [56] for a review). 
In general, these two effects are difficult to separate, making it necessary to model them 
jointly [4]. We thus include the effect of gravitational lensing in the shape correlations in 
our quantitative forecasts of Section 7. At linear order, the alignment (/) and lensing (G) 
contributions are simply additive, so that the full shape auto-power spectrum becomes 

C'yy(0 = / dkk^Pm{k)\Fl{k) + Fl^{k)\\ (4.20) 


where the lensing shear kernel for a distribution of source redshifts dNc/dz is given by (e.g., 
App. F of [51]) 


Ff^{k) = 


l{l + 2)\ 

2 (/ - 2) 


I rXraayi 

Uo 


dxk‘^ D^{k,z{x)) ^^^X 


f 




dz 


X 
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and is given by 


D^{k, z) = 




,,o (1 + z)D{z) 


A;2 


Z)(0) 


(4.22) 


The expression in Eq. (4.20) includes the intrinsic alignment auto-power spectrum, cross¬ 
correlations between weak lensing shear and intrinsic shapes, and cosmic shear. The cor¬ 
relation between intrinsic shapes and alignments occur when a galaxy in the background is 
being lensed by the same matter density field that tends to align the second galaxy [2, 57]. 
Similarly, the density-observed shape cross-correlation becomes 


= / k^dkPrn{k){Fl{k)^F^{k)\Fj^{k). (4.23) 

As this is an important point, we reiterate that the result in Eq. (4.20), with Eq. (4.14) 
and Eq. (4.21), is exact on linear scales, i.e. it does not receive further relativistic corrections.^ 
Eor the intrinsic contribution, we have argued this in Section 4.1. Eor the lensing contribution, 
Ref. [51] showed that the seemingly “Newtonian” expression for the shear in conformal- 
Newtonian gauge [Eq. (4.21)] is the proper relativistic expression for the observable shear 
at linear order. Note that this does not hold for the convergence k, where the observable 
magnification in fact receives several relativistic corrections. As discussed above, Eq. (4.23) 
assumes several simplifications in the number count kernel Ff'[k). 

Eor the numerical implementation, we use the Limber approximation for I > 60 for 
all angular power spectra, following [58]. At these multipoles, the Limber approximation 
provides an estimation of the angular power spectra with sufficient accuracy for our purposes, 
and yields better numerical convergence. We have verified that using the full integrals up to 

1 = 400 does not change our results. Eor the quadrupole {1 = 2), we use a value of k^m = 

2 X 10~® /iMpc”^ as minimum wavenumber, which matters for the shape auto-correlation as 
discussed above. Eor all other multipoles the value of /cmin is chosen small enough to ensure 
a converged result. 

5 Modeling the red galaxy population and alignments 

Red galaxies have been observed to be subject to intrinsic alignments with a scale dependence 
that is well approximated by the linear alignment model (Eq. 3.10) on scales above ~ 
10 Mpc/h [3, 8-10, 57]. On the other hand, late-type galaxies do not exhibit detectable 
alignment on large scales [11, 12]. Thus, in order to estimate the intrinsic alignment power 
spectrum as a function of redshift, it is necessary to model the red galaxy fraction. 

The fraction of red sources. 


/red(^) — 


dNj-ed/dz 

dNc/dz ’ 


(5.1) 


depends on the intrinsic properties of red galaxies (their spectral energy distribution and 
luminosity function) and on the properties of the survey (filter, magnitude limit, source 
redshift distrubution). We use the model constructed by Ref. [3] and applied by Ref. [59] 
to obtain the fraction of red galaxies and their mean i-band luminosity as a function of 
redshift for an LSST-like survey. The source redshift distribution expected for LSST, which 


'^This assumes that certain potential gauge artefacts have been removed from A 2 [40]. 


16 








Figure 4. The expected fraction of red galaxies observed by LSST as a function of redshift for a 
limiting magnitude of him = 25.3 (AB). The gray dashed curve shows dN^ed/dz with an arbitrary 
normalization. 

extends from z = 0.1 to 2 : = 3, is given in [60]; there are 26 galaxies per arcmin^ and 
we assume an i-band limiting magnitude of 25.3 (AB), consistent with the so called ‘gold 
sample’ of galaxies with signal-to-noise ratio > 20 and measured shears in LSST [61]. We 
assume a red galaxy luminosity function as observed by Ref. [62] from a combination of 
galaxy surveys between 0.2 < z < 1.2. In our model, the paucity of red galaxies at z > 1.5 
implies that intrinsic alignments have no contribution to the total shear above those redshifts. 
The redshift distribution of red galaxies in LSST using this model is shown in Figure 4. 
Throughout, we assume a survey with an area coverage of 18,000 square degrees, since the 
effect of masking is accounted for in the effective number density of galaxies. When projected 
over the range 2 = [0,1.5], LSST should observe 3 red galaxies per arcmin^. Notice that this 
is optimistic because we are not modeling any additional selection cuts on galaxy shapes that 
could potentially affect the red population. In our model, blue galaxies are assumed to have 
no alignments. 

We can relate the linear shape bias b{ to the convention commonly used in papers 
reporting on measurements of intrinsic alignment as follows: 

bl = —Aj[Lr)^rnClPcvitJJ^i (^-2) 

where 12^ is the matter density parameter today and Cipcrit = 0.0134 is hxed by convention 
to the SuperCOSMOS intrinsic alignment amplitude [63]. Later on, we will treat Cipcrit as a 
free parameter when making predictions of non-Gaussianity constraints from future surveys. 
Aj is a luminosity-dependent amplitude measured by Ref. [3] from a sample of Luminous 
Red Galaxies (LRGs) in the Sloan Digital Sky Survey [64], 

1 ^g + 0.25 

Ai{Lr) = (5.76;0:®0) , (5.3) 

where Lr is the average r-band luminosity at a given redshift and Lq is a pivot luminosity 
corresponding to an absolute magnitude of Mq(z) = —22 — 1.22: — 51og^Q h. We emphasize 
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that ongoing imaging surveys, such as the Kilo-Degree Survey^, the Dark Energy Survey^, 
the PAU survey^ and Hyper-Suprime Carr? will be able to test the assumptions on the 
dependence of 6{ with redshift and luminosity we have made here. 

Unlike in the case of the scale-dependent bias in galaxy number counts, ~ (6^ —l)5c, 
for which a fairly accurate estimate is known, no such estimate exists for the response of 
shapes to the power spectrum anisotropy quantified by b^Q. We will thus assume that 
is comparable to b{ at any given redshift, specifically 

^NG = ^NG ^ 

where we assume throughout = 1 but keep it explicit in order to illustrate the residual 
modeling uncertainty. 

Note that the effect of the luminosity dependence of the intrinsic alignment bias was 
not taken into account in Section 4.3 above, but will be included in the quantitative forecast. 
Its most important consequence is the existence of a Malmquist-like bias for the population 
of intrinsically aligned galaxies in a magnitude-limited survey. As a result, the red galaxy 
population at higher redshifts has a larger which results in a larger contribution from 
anisotropic non-Gaussianity compared to the case presented in Section 4.3, where b{ was set 
to a constant. 

6 Fisher matrix analysis 

In order to estimate the ability of future surveys to set joint constraints on isotropic and 
anisotropic non-Gaussianity, we perform a Fisher matrix analysis over the parameter set 
{6”, Aq, A 2 , Gipcrit, <78}- The remainder of the cosmological parameters are fixed to their 
values given at the end of Section 1. Further, we include a prior in cjs = 0.831 ± 0.013 from 
Planck [32, Table 3, column 4]. 

The predictions for the angular auto- and cross-power spectra were given in Section 4. 
All auto-power spectra also contain noise terms, which we assume to be white, 

= ^ and iV" = ^ , 
no rin 

where is the dispersion of the intrinsic shape including measurement noise per component, 
and fbG is the projected surface density of galaxies with shapes per steradian. Similarly, 
is the surface density of galaxies used to compute the clustering signal. 

Since all perturbations are Gaussian at the order we work in, the covariance between 
two power spectra is given by 

CoY[C^p{l),C,s{l')] = + [C^,{l)C08{l) + CMC0,{1)] ( 6 . 2 ) 

where the indices {a,/3,7, d} run over {re, G}, i.e. number densities and shapes. 

The Fisher matrix is given by [65, 66] 

®http://kids.strw.leidenuniv.nl 
®http://www.darkenergysurvey.org 

^“Physics of the Accelerating Universe”, http://www.pausurvey.org/home-PAU.html 
®http://www.naoj.org/Proj ects/HSC/ 


( 6 . 1 ) 


(5.4) 
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(6.3) 


F^ = Y. 

I 


dUil) 

dPf, 


Cov ^(/) 


dBjl) 

dpu 


where Cov is the covariance matrix of the data vector D = {Cnn-, Cn'y, and 5^ indicates 
the partial derivative with respect to the parameters of interest. We perform these derivatives 
analytically on the expressions for Ca/sil), and evaluate the results numerically. In practice, 
A 2 cannot be constrained independently of b^Q- While we assume b^Q = 1, this means that 
the power spectra are effectively constraining the product 5 ^qA 2 . The la uncertainty in 
each parameter, after marginalizing over all other parameters, are given by cr^^fun-margin = 
We also consider forecasted errors after marginalizing only over { 6 ”, Cipcrit; o's}; 
which allows us to determine whether degeneracies between ^0 and A 2 exist. If Aq and A 2 
correspond to coordinates I and 2 of the Fisher matrix, we obtain the partially-marginalized 
bounds from defining the following sub-matrix. 


{G ^)ij - [F 2]J=[1,2] 


and then the partially-marginalized uncertainties are given by 

^2, part—margin — 


(6.4) 


(6.5) 


7 Constraining non-Gaussianity with clustering and shapes 

In this section, we forecast the ability of galaxy imaging surveys to detect anisotropic non- 
Gaussianity using galaxy shapes, in particular those of red galaxies. We consider an LSST-like 
survey with 26 galaxies per arcmin^ [60] and we model the red galaxy population as described 
in Section 5. We assume that the dispersion in the intrinsic ellipticities per component is 
a.y = 0.3 [67]. We do not account for dispersion due to photometric redshift scatter in the red 
galaxy population in this section, since we are integrating over the full redshift range of the 
red galaxy population. Figure 5 shows the position-shape and shape-shape angular power 
spectra and their expected uncertainties as a function of multipole number for the cases with 
6 ^q ^2 = 0 and = 100. Clearly, the shape correlations are dominated by the leasing 

contribution, and the fractional effect of anisotropic non-Gaussianity is significantly smaller 
than in the intrinsic-only correlations shown in Figures 2 and 3. 

We combine the clustering, galaxy-lensing and cosmic shear information for the red 
galaxy population to obtain constraints on {Cipcritj (^ng^s); cs}- Notice that we 

assume that the luminosity and redshift dependence of the intrinsic alignment model is 
hxed by current observational data, but we allow the overall amplitude of alignments to be 
constrained through changing the normalization given by Cipcrit (no priors are considered). 
Analogously, we assume a fiducial constant, redshift-independent bias for red galaxies, 5” = 2 , 
that is allowed to vary in the Fisher analysis. We obtain the following fully marginalized la 
bounds on the uncertainties in the parameter set from the multipole range 2 < I < 600: 

A(Cipcrit) = 8.6 X 10 -^ 

A67 = 2.6 X 10“^ 

AAq = 44, and 

A( 6 ^gA 2 ) = 150. (7.1) 
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Figure 5. Position-shape correlations and shape-shape correlations of red galaxies expected to be 
observed in LSST between 0.1 < z < 1.5. The case of Gaussian initial conditons is shown in solid 
black, and the case with A 2 = 100 is shown as dashed gray lines. The error bars show the expected 
uncertainties for a Gaussian fiducial cosmology. 


Our forecast for the uncertainty in A 2 (assuming b^Q = 1) is somewhat worse than the current 
CMB constraints, AA^^^ ~ 90 [31, 68 , at 68 % confidence level]. We emphasize however that 
galaxy alignments probe different scales than the CMB, and thus these constraints should 
be seen as complementary to the CMB in view of a possible scale-dependence of the non- 
Gaussianity. This applies in particular to the small-scale modes, where the modes relevant for 
galaxy formation are ^5 > 1 /iMpc“^, while the CMB probes modes with ^5 < 0.1 /iMpc~^. 

Note further that the constraints depend on the fiducial value of A 2 chosen, where 
our choice of fiducial ^42 = 0 is conservative. This is because the terms oc (^ 2 )^ do not 
contribute to the constraints in this case. In particular, this applies to the quadrupole of 
the shape auto-power spectrum, which diverges as —)• 0 as discussed in Section 4.3, and 

is the most prominent non-Gaussian feature in Figure 5. Thus, the true log-likelihood is 
steeper than the quadratic approximation around ^2 = 0 suggests, so that even in case of a 
non-detection the actual constraints obtained for a survey with these specifications could be 
tighter. Further, in case of a detection of non-zero A 2 with shape correlations, the amplitude 
of C^^{1 = 2 ) could be used to constrain kmm via its logarithmic dependence on it. 

For the well-known isotropic non-Gaussianity on the other hand, our results indicate an 
uncertainty of = 11. This latter constraint is signihcantly worse than what is expected 

for Stage IV spectroscopic surveys, since we consider only a single redshift bin which loses 
information on all line of sight modes. We will return to this issue at the end of this section. 

Figure 6 shows likelihood contours in the ^0 ~ ^ng ^2 plane, obtained by marginalizing 
over Gipcrit, <78 and 6 ”. They correspond to { 1 , 2 , 3}(7 confidence intervals. We hnd that in¬ 
cluding all auto- and cross-correlations of number counts and shapes almost entirely removes 
any degeneracy between the two non-Gaussianity parameters. 

Given our linear treatment of the matter power spectrum, we consider the effect of 
restricting the forecasts to lower multipoles (note however that nonlinear corrections to the 
lensing contribution are straightforward to incorporate). Figure 7 shows the forecasted un¬ 
certainty in the anisotropic non-Gaussianity parameter as the maximum I is increased. This 
figure shows that although the effect of non-Gaussianity is most significant on large scales, 
the small scales still contain significant information in particular on 6 ” and Gipcrit which are 


- 20 - 







































300 



Figure 6. {1,2,3}(t uncertainties in Aq and b^QA 2 (thick to thin black, blue and red solid lines, 
respectively) marginalizing over Cipcrit, cg and 6”. These constraints are obtained from considering 
the combination of clustering, position-shape correlations and shape-shape correlations and using the 
redshift distribution of red galaxies expected to be observed in LSST between 0 < z < 1.5. 


helpful in improving the constraints on the non-Gaussianity parameters. While extending 
the modeling of the power spectra to nonlinear scales is beyond the scope of this work, there 
is thus potential for better constraints on the non-Ganssianity parameters from I > 600. 
However, the effects of baryonic physics on the nonlinear leasing and alignment power spec¬ 
trum on those scales are still a matter of debate (see, for example, [69]). On the other hand, 
removing scales below I < 10 while keeping Zmax = 600 fixed increases the uncertainty in 
the marginalized A( 6 ^qH 2 ) by a factor of 1.6. Not snrprisingly, the very low mnltipoles 
are important for constraining primordial non-Gaussianity, requiring careful control of the 
systematics in imaging surveys. 

In our constraints we have included the full redshift distribution of galaxies, thus sup¬ 
pressing any information contained in modes along the line of sight. This is the reason why 
our forecasted constraint on is ronghly an order of magnitude worse than what is obtained 
from spectroscopic surveys of similar volume. Given that, nnlike the leasing held, the tidal 
alignment is an intrinsically three-dimensional effect, we thus expect a potential signihcant 
improvement over the constraints quoted here by using a full tomographic analysis rather 
than a single redshift bin, which would restore information on the line-of-sight modes. This 
has been found to hold for the scale-dependent bias in the galaxy density indnced by 
[70, 71], and for lA in the case of Gaussian initial conditions [7, 72]. Further, tomography 
should help break any degeneracy between the alignment strength and cosmological param¬ 
eters such as (Is, which we have not considered here. In fact, a much smaller spectroscopic 
sample of red galaxies with well-measured shapes could deliver interesting constraints as well. 
We defer the forecast for the full tomographic and spectroscopic cases to future work. 

8 Discussion and conclusions 

We have shown that the intrinsic alignment of red (elliptical) galaxies will become an inter¬ 
esting probe of inflationary physics in the future. Through their sensitivity to a quadrupolar 
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Figure 7. Bound on A( 6 ^qA 2 ) after marginalizing over all other parameters and as a function of 
the maximum multipole considered (holding the minimum I fixed at 2 ). 


anisotropic mode coupling, galaxy shapes expand the parameter space of inflationary models 
that can be probed using large-scale structure observables. The anisotropic squeezed-limit 
bispectrum translates into a local anisotropy in the matter power spectrum, which imprints 
itself in the shapes of halos and galaxies forming at that location. This non-Gaussian bias 
produces a tendency towards radial or tangential alignments, depending on the sign of the 
non-Gaussianity, which increases rapidly towards large scales. This is the analog in galaxy 
shapes of the familiar scale-dependent bias in galaxy clustering induced by isotropic primor¬ 
dial non-Gaussianity. 

As a specific example, models such as solid inflation, the realization of inflation through 
the symmetries of a solid [22], could be probed through this method. More generally, as 
was recently pointed out by Ref. [26], the angular dependence of the bispectrum encodes 
information about the spin of light particles during inflation. Therefore, constraints on 
the quadrupolar part of the bispectrum, analyzed in this work, provide information about 
tensor degrees of freedom during inflation. Moreover, we have shown that should anisotropic 
non-Gaussianity be detected, the quadrupole of galaxy shapes can be used to constrain the 
maximum scale at which this mode coupling persists, which is presumably related to the 
beginning of or some other special event during inflation. This is a unique possibility of 
probing physics outside of our horizon, which is not possible with isotropic non-Gaussianity. 

We have modeled the red fraction of galaxies in the LSST survey following the approach 
suggested by [3], and conservatively assumed that late-type (blue) galaxies do not align, 
which is consistent with currently available constraints [11, 12]. We have shown that the 
constraints on A 2 for future surveys are comparable to, though slightly worse than, current 
GMB constraints. Specifically, we find AA 2 = 150(6^q)“^ for LSST, while current constraints 
on A 2 from Planck require —80 < A 2 < 110 at the 68 % confidence level. Note that this 
constraint was derived for a fiducial A 2 = 0, which neglects any terms scaling as (^ 2 )^. To 
achieve these constraints, imaging surveys will need to access scales below I < 10 , otherwise 
the constraint worsens by a factor of 1.6. However, a three-dimensional reconstruction of the 
auto and cross-spectra of positions and shapes could potentially yield significantly tighter 
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constraints on both A 2 and Aq, as the alignment process is intrinsically three-dimensional, 
in contrast to leasing. Just as importantly, the constraints on A 2 from galaxy shapes are to 
be considered complementary to those from the CMB, as they come from different physical 
scales {k > l/iMpc“^ \s k < 0.1/iMpc“^). 

One key caveat to these forecasts is the unknown value of the parameter which 
quantifies the response of galaxy shapes to an anisotropic initial power spectrum. In this 
work, we have assumed that is of the same order as the tidal alignment strength b\ which 
has been measured observationally. Ideally, should be estimated from hydrodynamic 
simulations of galaxy formation, by running a simulation with an anisotropic initial power 
spectrum and measuring the response of halos and galaxies to this preferred direction (also 
as a function of luminosity or mass). For the purpose of this work, we chose to simply phrase 
our constraints in terms of the combination where ~ 1 if our assumption is valid. 
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A Derivation of Eq. (3.11) 


The left-hand side of Eq. (3.11) is given in terms of the potential three-point function by 


(J(x)J(y)Ai,(y)) = I J^ J 


kiikij 1 
- 




(27r)3 ^ ^ J (27r)3 J (27r)3 kf 

X Al(/i:i)Al(A;2)(<('k0ki0k2), 

where r = x — y, and Ai is defined in Eq. (3.13). Using the definition of the bispectrum, 

{4’k<Pki<Pk2) = (27r)^j£,(k-F ki-F k2)R^(k,ki,k2), (A.l) 

we have 


(J(x)J(y)A:,,(y))= [ ^e^^-^Mik) ^ 


(27r)’ 


(27r)’ 


kl 


- -Si 


X M(ki)M(lk + kiDB^k, fei, |k + kil). 


(A.2) 


We now expand the integrand in powers of q = k/ki. Eurther, we define // = k • ki. As 
shown in App. C of [35], but generalized to the anisotropic case Eq. (1.1), the bispectrum is 
given by 

B^{k,ki, |k -h ki|) = '^AiVi{ij)P^{k)P^{ki) [2 (2g/r + [1 - 2f?]q‘^)n^] + 0{q^), 

i 


where = Ug — ^ and we have assumed a pure power-law P(j,{k) for simplicity. Eurther, 


M(|k + ki|) =M{ki) 


1 -h qnuMiki) + 0{q^) ; 


nM{ki) = 
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d\n.ki 


(A.3) 
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Inserting these expressions into Eq. (A.2), we obtain 


(S(x)S(y)Kij(y)} = 
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In the following we will neglect any 0{q^) corrections; in terms of observed correlations, these 
would lead to a small scale-independent correction which is degenerate with the Gaussian 
part. We now choose k to he along the z axis. We can then write 




^ COS^ (^(1 —//^) 
cos i^sin(/)(l — fj?) 
cos (f) qL\/l — /i^ 


cos (j) sin <^(1 — //^) cos 4> /iy^l — // 
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(A.4) 


Performing the (p integral only leaves the diagonal terms (as expected by symmetry): 


dp kiikij 

Jo 27r kl 


/1 - /i^ 0 0 \ 

0 I - /i2 0 

\ 0 0 2/i2 / 


(A.5) 


In order to obtain an expression for general directions of k, we notice that by symmetry, 
the final expression can only involve the tensors Pk^ and d^P Straightforward algebra then 
yields 
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We thus have 

{5{x)5{y)Kij{y)) = 
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In the last line, we have used the orthonormality of the Legendre polynomials. Note that 
since L is even, the qfi term does not contribute and corrections are order q^. The ki integral 
now simply gives the (formally divergent) variance of the density field ((5^). By dehning 

/ d^k 

-^e<^’^M(k)P4k) (A.7) 

we then obtain 

{6{x)6iy)Kij{y)) = ‘^A2Vij^^5{r){5‘^). (A.8) 

This is Eq. (3.11). It is easy to see that the corresponding result for the quadratic tidal term 
simply becomes 


(5(x) 


K,kK) - -6ij{Ki^y^ 


(y)) = ^A2VijC4>5{r){S'^) . (A.9) 
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This is because the analog of Eq. (A.2) for this term differs by factors constructed out of 
kij (k + ki)j/(/i;i|k + ki|), which at leading order are simply kukij/kf. 

Also, let us show that anisotropic non-Gaussianity does not lead to a scale-dependent 
bias in the number density of tracers. The tracer density locally depends on the variance 
of the density field smoothed on some scale, as well as the square of the tidal tensor Kfj 
[Eq. (3.2)]. Let us calculate these local quantities given the power spectrum Eq. (2.3): 


= I 


(27r)3 

(27r)^ 


1 -h /(ki • k 5 )(/)(kL)e*^^^ P{ks) 


tjg - 


1 -h /(ki • ks)4>{^L)P^^^ P{ks) 


(A.IO) 


Aligning with the z axis, we see that both terms involve the integral 



d^i fin) 


2Ao, 


(A.ll) 


i.e. these quantities pick out the monopole of /. Thus, non-Gaussianity of the quadrupolar 
type does not affect the local averages Eq. (A.IO) at leading order in the non-Gaussianity. 


B Renormalized shape biases 

In this section, we provide a rigorous derivation of the reasoning outlined in Section 3. Eor 
this purpose, we introduce a fictitious coarse-graining scale Rl which we assume to be much 
smaller than the scale on which we measure correlations. We denote coarse-grained fields by 
a subscript L. Since we will be dealing exclusively with galaxy shapes and not with number 
counts, we will drop the superscript I for the remainder of the appendix. 

Let us consider the case of Gaussian initial conditions first. Since Rl S> i?*, the shape 
of galaxies at position x becomes a local function of the coarse-grained density 5l and tidal 
field K^l ■ addition, the shapes are in general a functional of the small-scale fluctuations 
within a region of size ~ iZ* around x. Hence, 


g,^{^) = Ftj{6L{-^)-Kl\-K)-6s), (B.l) 

where 5s = d — 5 l are the small-scale fluctuations. The superscript L denotes the fact that 
depends on Rl- We can now perform a formal Taylor expansion of in 6l and K^l, 
taking into account that P/j has the same symmetries as gij. The first few terms of the 
expansion are 


5ij(x) = ciK]^{x.) + ic2(5L(x)A:];((x) -h 


Ki^K 


kj 


--^S'UK'Lf 


+ 0(«i, V^V). (B.2) 


The coefficients Cj depend on Rl and on the small-scale modes 5s- In the Gaussian case, 
the small-scale modes are uncorrelated with the long-wavelength perturbations, so that they 
merely lead to stochastic terms which do not correlate over large scales. Then, the Cj can 
be regarded as effective constants. We have not written the most general leading order term 
here, which would be CijkiK^\ and instead reduced Cijki to a scalar ci. The reason is that 
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in order to construct anything more general we would need some preferred direction which 
is absent in 3D space.® 

The expectation value of gij has to vanish by symmetry. Note that this is achieved 
at second order without having to add any counter term in Eq. (B.2). Now consider the 
transformation 


Kl{yi,T) K£{x,t) + 5l(x,t) 


(B.3) 


where D{t) is the matter growth factor, is constant and /3j* = 0. This corresponds to 
adding a uniform tidal field to the Newtonian potential, 


^'Ar(x) 


^'Ar(x) + + z)D{z)l3ijX’‘X^ . 


Under this transformation, the expectation value of g changes to 

{9ij{^,r))i3 = ci{T)(3ij + 0{Sl, V'^6 l) ■ 


(B.4) 
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In the following, we will assume we are evaluating the correlations at some fixed time t, and 
drop this argument for clarity. We now introduce the renormalized bias parameters through 




d{gij)fj 


dl3i. 


/ 3=0 


= ci + Oi[6L,K^if 


(B. 6 ) 


Note that the definition of bi is explicitly independent of Rl. The i?L-dependent quadratic 
corrections in Eq. (B. 6 ) come from cubic order terms in Eq. (B.2), hence we have not written 
them here. Similarly, the renormalized quadratic tidal bias is defined via 


bt^ 


d‘^{9ij)p 


diPijY 


/ 3=0 


= ct + Oi[6L,K^if 


(B.7) 


where now the quadratic corrections come from fourth-order terms in Eq. (B.2). 

In order to derive the renormalized bias corresponding to C 2 , we generalize Eq. (B.3) to 


(x,r) (x,t) +D{t) 13''^; + D{t) c, 


(B. 8 ) 


where the transformation of 6l is the same as the one used to define renormalized local bias 
parameters in [35]. Specifically, it corresponds to adding a uniform matter density component 
cp. We then define 


b2 


d‘^{9ij)p. 


dl3ijdc 


13=0,c=0 


C2+0{[6L,K^i?). 
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Again, 62 is independent of i?L by construction. 

We now turn to correlations of gij on large scales. Eor simplicity, let us consider the 
cross-correlation with matter here, i.e. {6{x)gij{y)). After renormalizing the Cn into bn we 
obtain 


( 5 (x) 5 *j(y)) = 6i(5(x)A:];((y)) -62 <^<5(x)<5L(y)A:]f (y) 


+ lbt 


<5(x) 
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®In reality, the line of sight n provides such a preferred direction. Thus, in principle a linear term like 
n''h^hkfut'^ could be present due to redshift space effects. However such a term will not be relevant for 
photometric galaxy catalogues. 
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where the ellipsis stands for higher order no-zero-lag correlators, i.e. correlators that do not 
contain any zero-lag terms such as {KijK^^) when expanded into connected correlators. 
If |x — y| 3> Rl (and the correlation function does not have sharp features), then 

(y)) = =Pij^(|x-y|) (B.ll) 

where ^{r) is the matter correlation function. This renormalized bias quantifies the linear 
response of galaxy shapes to a long-wavelength tidal field as described above. For Gaussian 
initial conditions, the three-point function in Eq. (B.IO) is only due to non-linear evolution 
and relevant on smaller scales. 


B.l Non-Gaussian initial conditions 

We now turn to the case of non-Gaussian initial conditions, quantified by Eq. (1.1). Note 
that the three-point term in Eq. (B.IO) is sensitive to the bispectrum in the squeezed limit if 
“ y| Rl, the coarse-graining scale adopted. As shown in App. A, for a squeezed-limit 
bispectrum of the form Eq. (1.1) we have 

{s{:x.)6L{y)K^i (y)) = ‘^A 2 VijCs 4 >{\^ - y|)(^i, (B.12) 

where is the matter-potential cross-correlation and (t£ is the variance of the coarse¬ 

grained density field, a‘j^ = (Sf). A similar result (with prefactor 2/15 instead of 2/5) holds 
for the term oc bt- The fact that this term scales as i.e. that it strongly depends on the 
arbitrary coarse-graining scale indicates that the bias expansion Eq. (B.2) is not sufficient in 
the presence of primordial non-Gaussianity. This exactly parallels the case of galaxy density 
correlations, where the corresponding three-point term also scales as a‘j^ [35]. 

In order to remedy this unsatisfying situation, we explicitly include the dependence of 
Qij on the small scale fluctuations indicated in Eq. (B.l). Of course, this dependence has 
to be consistent with the symmetries of gij. Physically, we expect that the lowest order 
dependence on 5s in Eq. (B.l) will be through a quadrupolar asymmetry in the statistics 
of the local small-scale density field; an isotropic change cannot produce an effect at linear 
order because we are expanding a spin-2 function. Note that it is sufficient to phrase the 
dependence on the small-scale modes in terms of 5s, since the small-scale tidal field is related 
to 5s through the Poisson equation. Let us define the parameter as the anisotropy of the 
local small-scale correlation function within a region of size R^: 


ys''{^) = ^ [ d^rWL{\r\)5six)T>'^''5six + r) 

/ ySjL 

-^WL{ki)W^{ki)P^{k ^), (B.13) 


where VFi is an arbitrary isotropic filter function on the scale Rl (e.g., Gaussian or tophat), 
and Wl is its Fourier transform. denotes the derivative operator Eq. (3.3) acting with 
respect to r. While this definition is convenient, other definitions are certainly possible and 
do not change the results. Note that since by definition 5s = 5 — 5 l, the statistics of 5s also 
depend on Rl- We now let gij depend locally on 

yij(x) = Flj{5L{^)-, ATi (x); yf (x)). (B.14) 
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We expand Fij to linear order in which adds a term 


CNGVs (x) 


(B.15) 


to the expansion on the r.h.s. of Eq. (B.2). At linear order in the non-Gaussianity parameters 
Ai, we do not need to consider higher order terms involving ylF The same symmetry argu¬ 
ment as for dictates a simple scalar coefficient cng- Correspondingly, we have another 
contribution to the matter-shape correlation: 




(y)) + 
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+ CNG((^(x)y*^(y)) H- 


(B.16) 


The last term can be calculated along the same lines as the calculation in App. A. The 
Fourier transform of Eq. (B.13) is given by 
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We then obtain 
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where Ais{k) = Ai{k)Ws{k). The derivation now exactly parallels that of Eq. (A.8), where 
the ki integral now yields 


/ ^^WL{ki)Ml{ki)P^{ki) = , 

cancelling the prefactor. We obtain 

{S{y)yl^{^)) = ‘^A2Vij^^s{r) ■ 


(B.19) 


(B.20) 


There are thus three non-Gaussian terms that are proportional to A 2 'Dij^Scj)^ two of which 
are strongly dependent on Rl. Intuition tells us that these contributions should combine 
to an i^L-independent term involving a renormalized bias parameter 6ng- As before, 6ng 
should be given as the response of the mean shape of galaxies to a specific (iiL-independent) 
transformation of the density field. Consider the transformation 
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or equivalently in Fourier space 
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This corresponds to an anisotropic, scale-independent rescaling of the density field. In the 
following we will assume without loss of generality that aim is trace-free. The expectation 
value of the shape after this transformation is 


(9«>« = \l>2{SLK‘i)c + it. ^ + cng(!/?). 


(B.23) 
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After some algebra, this yields 
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The linear response of the mean shape of galaxies, which shall define our bias 6 nG) is thus 
given by 


f^NG = 

Solving this for cng yields 
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which we insert into Eq. (B.16) to obtain 


{5{-x)gij{y)) = hiVij£,ss{\^ - y|) + 36NG^2Aj?5<A(|x - y|). (B.28) 

Our intuition is thus confirmed: by adding the bias 6 ng, which quantifies the linear response 
of galaxy shapes to an anisotropic power spectrum of initial fluctuations, we have absorbed 
the unphysical dependence on a\ present in the local quadratic bias expression. 


C Non-Gaussian contributions to weak lensing shear 

In Section 4 and following, we have neglected any modification of the lensing contribution to 
galaxy shapes due to primordial non-Gaussianity. We now show that there is a modification, 
but it is negligible compared to the scale-dependent signature in the intrinsic alignments. 
The observed shear is not exactly given by its linear expression 7 introduced in Section 4, 
but rather there are nonlinear corrections, 

7°bs(0)=^(0)+C.,N(0)7(0) + -- - , (C.l) 

where is a constant of order unity. The nonlinear corrections are due to reduced shear [56] 
as well as lensing bias corrections induced by source density weighting [74]. This leads to a 
three-point contribution to shear two-point functions of the form 

(^obs(^)^obs(0/)) ^ (^(0)^(0/)) + 2c^{yeueue')) + • • •. (c.2) 
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Primordial non-Gaussianity in general contributes to this three-point function. This was 
studied in detail by Ref. [47], who showed that 


/ j2; 

(27r)2 |1 - li|,Z). 


(C.3) 


Here, is the convergence bispectrum, and we have aligned li so that (fn = 0. Further, we 
have used that in Fourier space 7(1) = 71 ( 1 ) -|- *72(1) = exp(2i(/);)«;(l). Assuming the Limber 
approximation, the convergence bispectrum is related to the matter bispectrum by 
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where VV’l(Xs, x) = x/XsiXs—x) denotes the lensing kernel. Ref. [47] only considered isotropic 
local non-Gaussianity, however their results are easily generalized to the anisotropic case by 
making use of the squeezed limit. While in the isotropic case, the low-/ contribution to 
Eq. (G.3) is canceled by the cos2(/>7 modulation, this factor will lead to a non-zero contribu¬ 
tion in the anisotropic case, which will be of the same order as the isotropic contribution to 
C. obs (/) discussed in [47]. 

For our purposes, an order of magnitude estimate suffices. By using the squeezed limit 
for the bispectrum and approximating the projection kernels, we can derive in analogy with 
Sec. II in [47] 
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Thus, at the values for A 2 relevant for detectable alignment signatures, the fractional modi¬ 
fication to the shear is much less than one even at the lowest multipoles. 
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